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Field and laboratory measurements of suspended sediments over wave ripples show, for time- 
averaged concentration profiles in semi-log plots, a contrast between upward convex profiles for 
fine sand and upward concave profiles for coarse sand. Careful examination of experimental data 
for coarse sand shows a near-bed upward convex profile beneath the main upward concave profile. 
Available models fail to predict these two profiles for coarse sediments. The 1-DV gradient diffusion 
model predicts the main upward concave profile for coarse sediments thanks to a suitable /3(y)- 
function (where /3 is the inverse of the turbulent Schmidt number and y is the distance from the 
bed). In order to predict the near-bed upward convex profile, an additional parameter a is needeed. 
This parameter could be related to settling velocity (a equal to inverse of dimensionless settling 
velocity) or to convective sediment entrainment process. The profiles are interpreted by a relation 
between second derivative of the logarithm of concentration and derivative of the product between 
sediment diffusivity and a. 



I. INTRODUCTION 

The description of sediment transport in the marine environment is very complex due to complex interactions between waves, 
currents, sediments and sea-bed. The understanding of the physical processes of sediment dynamics in this environment is of 
crucial importance for accurate predictions of near-shore sediment transport and sea bed changes. In pure wave or oscillatory 
flows, turbulence is restricted to a thin boundary layer. The presence of a current together with waves implies that the 
turbulence spreads to cover all the flow depth. Very close to the bottom, the turbulence is due to the oscillatory motion and 
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predicted concentrations do not differ very much whether the turbulence due to the current is included or not 9]. Suspended 
sediment concentrations depend on the interaction between suspended sediments and the flow's turbulence. The feedback 
between sediment and turbulence has an important effect on sediment transport 5|. When considering suspended sediment in 
oscillatory flows, it is necessary to distinguish between the sheet flow regime (Shields numbers > 0.7) and the ripple regime at 
lower shear stresses. In extreme wave conditions and/or shallow water the sea bottom is plane and sheet flow is the dominant 
transport mode !§]]. However, for moderate wave conditions and/or deep water, wave ripples can be formed on the sea bottom. 
If the ripples are relatively steep (r) r /\ r > 0.12, where r) T is the ripple height and A r is the ripple wavelength), the mixing 
close to the bed is dominated by coherent, periodic vortex structures 14] 3] 1^| 18]. Above rippled beds, the mixing in the 
near-bed layer is dominated by the mechanism of vortex shedding which entrains sediments near flow reversal during each wave 
half cycle to several ripple heights 2(| . 

Examination of concentration profiles for different sand sizes suspended in the same flow over wave ripples provides infor- 
mation about the involved mechanisms. Field and laboratory measurements of suspended sediments over wave ripples show a 
contrast between an upward convex concentration profiles, time-averaged in semi-log plots, for fine sand and an upward concave 



profiles for coarse sand 
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The classical gradient (Fickian) diffusion model allows a good description of fine sediment 

n n 

concentrations but fails to describe concentration profiles for coarse sand [la]. The finite- mixing- length model [la], which is of 
similar generality to the Lagrangian Taylor's model, contains gradient (Fickian) diffusion model as a limiting case for l m /L — > 0, 
where l m is the mixing length and L the scale of the distribution under consideration. Even if the finite-mixing-length model 
improves the concentration profile c(y) for coarse sand (where y is the distance from the bed), a difference with experimental 
data remains. However, it is possible to improve this profile by a suitable turbulent mixing which is able to increase the upward 
concavity of c(y) [lj. This is possible with a turbulent diffusivity e s which takes into account the effect of grain size. The 
diffusivity of sediments e s is related to the diffusivity of momentum, the eddy viscosity ut, by a coefficient j3 — e s /v t (i.e., the 
inverse of the turbulent Schmidt number). In order to include the effect of sand grain size on e s , it is important to consider a 
realistic function for /3 which is dependent on grain size and concentrations and therefore on the distance from the bed y. 

The value of /3 has been the subject of much research. Researchers found that p approaches unity for fine sediments and 
deviates for coarse ones 8] ly^. In suspension flows over movable beds, experiments show that depth-averaged /3- values are 
smaller than unity (/? < 1) for flows without bed forms while they are larger than unity (/? > 1) for flows with bed forms 1C|. 
In oscillatory flows, it is known that cycle-mean sediment diffusivity above ripples is significantly greater than the cycle-mean 
1 3 18]. 



eddy viscosity, i.e. f3 > 



181 ]. This difference has not yet been clearly explained, but seems rest on the spatial-temporal 



correlation between suspended concentrations and vertical velocities in the flow field [18j. For 1-DV models, the value of /3 
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However, /3 should depend on grain size and 



was suggested empirically a constant equal to about 4 for rippled beds 
concentrations and should be therefore y-dependent. 

Careful examination of measured concentrations for coarse sand shows a near-bed upward convex profile beneath the main 
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upward concave profile. This near-bed profile was related to particles settling velocity which decreases near the bottom for 
high concentrations. Experiments have demonstrated that particle settling velocities are lower at higher concentrations. This 
behavior is described by the well known Richardson and Zaki 17| equation which depends on an empirically determined 
exponent. This equation was adapted for natural marine sands by Baldock et al. 2]. 

The aim of this study is to: 1) improve the prediction of concentration profiles by the gradient (Fickian) diffusion model; 2) 
provide an analytical interpretation of the main upward concave concentration profile and the near-bed upward convex profile 
for coarse sand. This interpretation requires a relation between sediment diffusivity, particles settling velocity and sediment 
concentrations. The following Section presents gradient (Fickian) diffusion and finite-mixing-length modeling of suspended 
sediment concentrations. The effect of sand grain size on sediment diffusivity and concentration profiles will be considered in 
Section 3. Section 4 presents possible effects of settling velocity or convective sediment entrainment process on the near-bed 
upward convex profile for coarse sand. 



II. GRADIENT (FICKIAN) DIFFUSION AND FINITE-MIXING-LENGTH MODELING OF 

SUSPENDED SEDIMENT CONCENTRATIONS 

Time-averaged concentrations of suspended sediment result from the balance between an upward mixing flux g m and a 
downward settling flux q" s = c(y)uj s as 

9m - c(y)ui a = (1) 

where to s is the particle settling velocity. 

A. Classical gradient (Fickian) diffusion model 

The gradient or Fickian diffusion model assumes that the mixing flux is proportional to the concentration gradient 



q™ = -Es-r- (2) 
dy 

where e s is the sediment diffusivity. Equations ^ and ((2)1 allow to write the classical diffusion equation 



d c 
dy 



(3) 
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Equation ([3| needs the sediment diffusivity e s which is related to the eddy viscosity vt by the parameter ft (i.e., the inverse 
of the turbulent Schmidt number) 

e s = P v t (4) 

In this equation, /? describes the difference between diffusivity of momentum (diffusion of a fluid "particle") and diffusivity of 
sediment particles. 



In the finite-mixing-length model 



B. Finite-mixing-length model 

lfj the swapping of fluid parcels (including suspended sediment) between different levels 
generates a net vertical flux of momentum and suspended sediment. The sediment concentrations in a lower parcel and in an 
upper parcel are respectively c(y — lm/2) and c(y + l m /2). If the parcels travel vertically with equal and opposite velocities 
±u m , the resulting flux is 

q m = v m [c(y - l m /2) - c(y + l m /2)] (5) 
where v m is the mixing velocity and l m the mixing length. The Taylor expansion of c(y — lm/2) and c(y + l m /2) gives 



v m l r , 



dy + 24dy 3 



Inserting Eq. ((6]) into |T]) and by including only the first two terms of the Taylor expansion in the brackets 



order ordinary differential equation (for more details see [lg|) 



3) 



3 



(6) 

gives a third 



l m d c dc uj s 

24 Ty~5 + + ZJZ C(y ' = 



(7) 



Equations (J2J and ([6]) give 



1 + 31 



d 3 c 
l m dy 3 



24 dc 
dy 



+ 



(8) 



If the eddy viscosity v t is given by a similar expression to ([8| with c replaced by u, the parameter /3 is determined by the 
ratio between the function in the brackets for each quantity or with L~ 2 ~ (d 3 c/ dy 3 ) / {dc/ dy) and L~ 2 w (d 3 u/dy 3 )/(du/dy) 
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8 * 1+l f^ + 

V t , 'm r— 2 



(9) 



C. Turbulent mixing 



In this study we will use a simple algebraic eddy viscosity model given by 



-Ci 



y 



ft = ai ft u„ y e 



(10) 



where k is the karman constant (= 0.41), it* the friction velocity, 5 the boundary layer thickness and ai and C\ two parameters 



n 



Eq. (|10[) was used in oscillatory boundary layers over rough flat beds by Hsu and Jan ll[ and over wave ripples by Nielsen 
and Teakle 

Hsu and Jan [llj showed that the vertical profile of Eq. (|10|l is similar to the time-averaged k-e two-equation model. Predicted 
velocity profiles in oscillatory boundary layers showed very good agreement with experimental data 11 1. 
Nielsen and Teakle 16(] used a similar eddy viscosity formulation given by 

{y - yo) 



v t = Ay v m o e - 



Lv 



(11) 



where v m o = v m (y = yo) is the mixing velocity at level yo; X and L v are two parameters. This model is based on a linearly 
increasing mixing length and an exponentially decreasing mixing velocity given respectively by 



lm — Xy and v m = v m0 e 



(y - yo) 



(12) 



The exponentially decreasing mixing velocity v m was chosen in agreement with turbulence measurements of du Toit and 
Sleath _ 



?} who observed that the turbulence intensity decayed exponentially upwards from a level yo very close to the ripple 

yo ' 

and L v 



crest. Eq. (ITTjl reverts to Eq. (fTO)) with ai k u* = A v m o e 



D. Gradient (Fickian) diffusion versus finite-mixing-length model 



In order to analyze concentration profiles obtained by gradient (Fickian) diffusion and finite-mixing-length models, we use 
experimental data obtained by McFetridge and Nielsen 
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1. Experiments of McFetridge and Nielsen 



In these experiments 



Q, 



suspended sediments are due to non-breaking waves over rippled beds. In order to establish reliable 
time-averaged concentration profiles, an extensive series of trials under a single set of wave conditions was chosen rather than 
a limited number of trials under a variety of wave conditions. Natural beach sand was used, wave maker piston amplitude, 
period, and water depth were adjusted to maximize near bed velocities and therefore suspended sediment concentrations. 
The flow parameters are: wave period T — 1.51 seconds, mean wave height H — 13 cm, maximum near-bed flow velocity 
Uom = 27.8 cm/s, near-bed flow semi-excursion a — 6.68 cm and maximum return velocity at the bed = 21.6 cm/s, mean depth 
of flow h — 30 cm. The ripples which were produced were highly uniform and regular with a mean ripple height of 1.1 cm, a 
mean ripple length of 7.8 cm and therefore a mean ripple steepness of 0.14. Measured concentrations were obtained by sieving 



;ccpn 

suction samples from different elevations above ripple crest 131 ]. 



2. Time-averaged concentration profiles 

Fig. (1) shows a comparison between time-averaged concentration profiles, ln(c) versus linear y plots, obtained with gradient 
(Fickian) diffusion and finite-mixing-length models and measurements ]j| for fine (d = 0.06 — 0.11 mm) and coarse sediments 
(d = 0.42 — 0.60 mm), where d is grain diameter. Settling velocities are respectively uj s — 0.65 cm/s for fine and ui s — 6.1 cm/s 
for coarse sediments 

Predicted concentration profiles were obtained by solving differential equations <[3j and (0 with respectively (|lip and (I12|l . 
We use the parameters proposed by Nielsen and Teakle 16j, namely A = 1, yp = 0.005 m, v m o = 0.025 m/s, L v = 0.022 m. 
Concentration profiles (figure 1) are the same as profiles of figure (6 of 16j). Even if the upward convex profile for fine 
sand versus upward concave for coarse sand is reproduced with the finite-mixing-length model, it remains a difference with 
experimental data for coarse sand at the top of Fig. (1) (y between 0.04 and 0.09m) which seems to become more important for 
y > 0.09 m. Since the aim of the study is to predict observed differences between concentration profiles for fine versus coarse 
sand suspended in the same flow, the turbulent flow's parameters were the same for fine and coarse sand (figure 1). 

III. EFFECT OF SEDIMENT DIFFUSIVITY ON CONCENTRATION PROFILES 

In order to investigate the effect of sediment diffusivity on concentration profiles, we consider in this section a settling velocity 
constant over the boundary layer thickness, equal to uj s t the terminal settling velocity in an infinite fluid. In our study we will 
consider that (3 approaches unity for fine sediments and deviates for coarse ones. 
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o a 



c/c Q [vol/vol] 



FIG. 1: Time-averaged concentration profiles over wave ripples: gradient (Fickian) diffusion model versus finite-mixing-length 
model, symbols, measurements 13j |. o, Fine sand (cj a = 0.65 cm/s); x, coarse sand (oj s = 6.1 cm/ a). Dash-dotted lines are 
solutions of gradient (Fickian) diffusion model (Eq. [3| while solid lines are solutions of finite-mixing-length model (Eq. 0, 
A = 1, y = 0.005 m, v m0 = 0.025 m/s, L v = 0.022 m. 



A. Concentration profile for fine sediments 

The assumption ft « 1 for fine sediments allows us to write 



i/i = aiKU,je 



-Ci 



v 



(13) 



Careful examination of time-averaged concentration profile (Fig. f) for fine sand {ui s = 0.65 cm/s) shows a small near-bed 

upward concave profile followed by the main upward convex profile (Fig. 2. a). These profiles are associated to a sediment 

diffusivity e s profile which increases then decreases (Fig. 2.b). The decreasing sediment diffusivity e s allows an upward convex 

concentration c(y) profile while an increasing sediment diffusivity e s allows an upward concave c(y) profile. 

d 2 In c 

Upward concavity /convexity of concentration profiles is, on ln(c) versus linear y plots, related to — — — , while increas- 
ed y 

ing/decreasing of e s is related to —r—. In order to understund this link, we need therefore a relation between — ; — ^ and 

dy dy 2 

de 3 



dy 
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FIG. 2: Upward concavity /convexity of concetration profile for fine sand versus increasing/decreasing sediment difusivity. (a) 

Upward concavity /convexity of c(y) (related to — ; — — ). (b) increasing/decreasing e a (related to — ■). 

dy dy 

We are able to write 

din c f / dc\ , . 

-W = -c\j-y) (14) 

and with Eq. fl3) 

dine uj 
dy e. 

The derivative of Eq. ([3)l and Eq. (|15|l are respectively 

d 2 c u) a c ( d e. 

and 



(15) 



dy 2 el \dy +U ^ (16) 



1 


' d 2 c 


1 


d c 


c 


dy 2 _ 


_ d 2 


dy 



(17) 



B Concentration profile for cdilse ESEMQiEsOF SEDIMENT DIFFUSIVITY ON CONCENTRATION PROFILES 
Inserting Eq. © and Eq. (fT6|) into Eq. (fTTj) . we obtain 



(18) 



dy 2 



1 


UJ S C 


( ds a 




1 




c 




\dy ' 


f U S ^J 


~ ? 





After simplifying Eq. (|18|) . we obtain 



d Inc w» d e a 



dy 2 e 2 a dy 



(19) 



Eq. (1191) provides a link between - } nC and — and therefore between upward concavity /convexity of concentration 

dy dy 

profiles and increasing/decreasing of sediment diffusivity profiles. Since ^ is always > 0, - ^ U C and 'LEl have the same sign 

si dy 2 dy 

and therefore increasing sediment diffusivity allows upward concave concentration profile, while decreasing sediment diffusivity 
allows an upward convex concentration profile (Fig. 2). 



B. Concentration profile for coarse sediments 



Due to grain size of coarse sediments, e s becomes different from v t . The difference in the shape of concentration profiles of 
fine and coarse sediments seems to be related to sand grain size which impacts sediment diffusivity through the parameter /?. 

The analysis presented in the last section suggests that an increasing sediment diffusivity is needed to predict the upward 
concave concentration profile for coarse sand. 

We write a /3(i/)-function as 



P = faU{y) 

where fa is the value of /3 close to the bed. 

In order to obtain an increasing sediment diffusivity, we suggest 

Cg- 

fe(y) = e s 



(20) 



(21) 



where Cp is a coefficient which should depend on sand grain size. 

No experimental data of /3(y) above wave orbital ripples are available to validate the proposed /3-function. Available data 
from measurements, in steady and uniform open channel flows with Acoustic Particle Flux Profiler lOj, show that ft increases 
with the distance from the bed, up to a given distance, indicating the higher the concentration, the lower are the /?- values. The 
proposed /3-function seems to be in agreement with these experimental data (fig. 3). 

It is possible to obtain a theroretical f3(y) profile, based on finite-mixing-length model, from Eq. ((9} since Z m is y-dependent. 
With a linear mixing length l m = A y, Eq. © gives 



P = 1 + C P2 y 2 



(22) 
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B Concentration profile for coarse sifflmfflBE NEAR-BED UPWARD CONVEX PROFILE FOR COARSE SAND 




FIG. 3: Vertical distribution of dimensionless coefficient /3/f3b- o, experimental data 10] (flow-sand condition Q70S'025II). 
Curves, solid line: Eq. (|21|) . dashed line: (Eq. 1221) from finite-mixing-length model with l m = Ay. 

where C732 is a coefficient (C$i = (A 2 /24) (Lj 2 — L^ 2 )). Fig. (3) shows that the profile of Eq. (|22[) (dashed line) confirms 
the shape of Eq. (|21|) (solid line) and suggests that fti, = 1. 

Using the ,3(i/)-function (Eq. 121 \ and eddy viscosity (Eq. I10|) . the sediment diffusivity is therefore given by 



fa ^ -(Ci -^)f 



a a y e 3s 



(23) 



where a 3 — fit, cti k u, and 5 a 



(Ci - Cp) 



are the two parameters of sediment diffusivity. 



Results of time-averaged concentration profiles over wave ripples are presented in Fig. (4) for fine {uj s — 0.65cm/s) and coarse 
(lo s — Q.lcm/s) sand. Comparisons with experimental data 13J show good agreement. The upward convex profile for fine sand 
versus upward concave for the coarse sand is reproduced with the gradient (Fickian) diffusion model. The imperfection which 
was observed in figure (1) for coarse sand (y between 0.04 and 0.09 m) is improved. 

Our results show that the decreasing sediment diffusivity e s for fine sand allows an upward convex concentration c(y) profile 
while an increasing sediment diffusivity e s for coarse sand allows an upward concave c(y) profile (Fig. 4). 



IV. THE NEAR-BED UPWARD CONVEX PROFILE FOR COARSE SAND 



Careful examination of measured concentrations for coarse sand, in semi-log plots, shows a small near-bed upward convex 
concentration profile beneath the main upward concave profile (Fig. 5. a). Examination of other experimental data 8] confirms 
the presence of near-bed upward convex profiles (Fig. 5.b). 
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c/c [vol/vol] 



FIG. 4: Time-averaged concentration profiles over wave ripples: gradient (Fickian) diffusion model with proposed /3(j/)-function 



symbols, measurements 



Fine sand (uj s = 0.65cm/s); x, coarse sand (uj s = Q.lcm/s). Curves, model with a a — 0.025m/s; 

= 0.24 m. 



S S 

dashed line (fine sand) S s = L v — — = 0.022 m; solid line (coarse sand) 5 3 — — — — - 

Ci (Ci — Cp) 



The 1-DV gradient diffusion model predicts the main upward concave profile for coarse sediments thanks to an adequate 
/3(j/)-function. The resulting sediment diffusivity e s with a constant settling velocity (w s = u s t) is unable to predict the near-bed 
upward convex profile. In order to predict this profile, we need an additional parameter. The sediment diffusivity should be 
modulated by an ^-dependent parameter a which depends on sand grain size and decreases with y. If we consider a constant 
settling velocity (uj a = oj s t), we write equation ([3]) as 



— = — c (24) 

dy ot(y)e a (y) 



Eq. (|19[) becomes therefore 

d 2 Inc _ lo s d (a e s ) 



(25) 



dy 2 (ae s ) 2 dy 

As in the analysis presented in the last section, Eq. (|25|) provides a link between upward concavity /convexity of concentration 
profiles and increasing/decreasing of a s s . Increasing a e s allows upward concave concentration profile, while decreasing 
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FIG. 5: Examination of experimental data; near-bed upward convex concentration profile, (a) Upward concave and near-bed 
upward convex concentration profiles (solid line) for coarse sand (data [13j]). (b) Time-averaged concentration profiles (data 
); o, Fine sand d=0,13 mm ; x, coarse sand d=0,32 mm; near-bed upward convex profiles (solid lines). 



a e s allows an upward convex concentration profile. Since e s increases with y for coarse sand, the near-bed upward convex 
concentration profile needs a decreasing a. 

A. Possible effect of settling velocity 

The near-bed upward convex profile could be related to settling velocity which decreases for high concentrations. Experiments 
have demonstrated that particle settling velocities are lower at higher concentrations. This behavior is given by the well known 
semi-empirical Richardson and Zaki I7| equation. We write Richardson and Zaki's equation as 



(26) 



where Co = c(y = yo) and n is an empirically determined exponent dependent on the particle Reynolds number R t at u s t and 
is constant for a particular particle. This exponent was determined experimentally as between 4.65 and 2.4 for increasing Rt. 
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A review of empirical expressions for n is given by Di Felice [6(]. These empirical values describing the variation of n with 
Rt were derived from fluidization and sedimentation experiments with spheres. Cheng [4J] proposed expressions for n which 
are functions of both Reynolds number and concentration. A method of determining n for natural sands was proposed by 
Baldock et al. [2| and was verified against laboratory fluidization data for marine sands. Settling velocity could be therefore 
not a constant on the vertical, and could depend on sediment grain size and concentrations and therefore on y. We write an 
{/-dependent settling velocity, required by our analytical study, as 

<jj s (y) = w a t f a (y) (27) 

where f s (y) is a function which is equal to 1 far from the bed where concentrations are small (c fn 0) and decreases near 
the bottom for high concentrations. In Section 3, we assumed that ui s = uj s t = constant and therefore f s = 1. However, if 
this approximation is valid only far from the bed where concentrations are very small, the observed near-bed upward convex 
concentration profile for coarse sand could be related to f 3 (y)- We need a suitable function for f s (y) in order to evaluate the 
effect on concentration profiles. 

1. y -dependent functions for settling velocity 
We propose two {/-dependent functions for cj s . The first is based on a simple concentration profile (Fig. 6. a) given by 



co V L, 

Inserting equation (|28|) into equation (|26[) . we obtain an {/-dependent function for settling velocity 
The second is an empirical function given by 

fs{y) = ^ = ^ 7 ^ 7 - (30) 

lis 

where y s and h s are two parameters which should depend on concentrations and particles size. In order to validate these 
functions, data of settling velocity w s were derived from measured concentrations for coarse sand and Richardson and Zaki's 
equation. Data of f s (y) are therefore obtained from equation (|26[) and measured c/co values, as 

«»>—(^) =M^) )" (3i) 

\ ' exp \ \ u / exp I 

Experimental data obtained from equation (|31[) depend on the value of the exponent n. Fig. (6.b) shows that the second 
function seems to be more accurate. 
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FIG. 6: (a) Concentration profile for coarse sand given by Eq. (I28[l . (L c ~ 0.008 m); X: measurements, (b) Settling velocity 
functions; x : measurements from Eq. (|31|l with n = 3.4; dashed line: first settling velocity function (Eq. 1291 L c = 0.008 m, 
yo = 0.005 m and n = 3.4); solid line: second settling velocity function (Eq. [30l y s = 0.018 m and h s = 0.0035 m). 

2. Effect on concentration profiles 

Fig. (7) allows to understand the possible effect of f s (y) on concentration profiles. In Section 3, f 3 was taken to be equal 

to 1 and concentration profiles depend only on s 3 . Here, f s is equal to 1 only far from the bottom where the concentration is 

small (Fig. 6.c). Fig. (7.a) shows l/f s and Fig. (7.b) e s /f B . The shape of e 3 /f s allows to predict the near-bed upward convex 
concentration profile for coarse sand (Fig. 8). 

Taking into account (Eq. 1271) . equation ([3]) becomes 

dy e s {y) s s (y) e s {y) 

My) 

In Eq. (|32[) . a — l/f s . Concentration profiles depend therefore on the ratio between sediment diffusivity and dimensionless 
settling velocity ^, s |^j and not on the sediment diffusivity e s (y) alone. 
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FIG. 7: Effect of dimensionless settling velocity profile fs(y) on predicted concentrations, (a) l// s . (b) e s /fs- 
Eq. (|19[) becomes therefore 

a 1 ] m r, ,., . d \ T ) 

(33) 



dy 2 fe^ 2 

U 

We are able to generalize the result found in Section 3 by: upward concavity/convexity of concentration profiles is, on ln(c) 

£ fy\ d 2 In c f (y\ 

versus linear y plots, related to increasing/decreasing of , . Eq. (1331) provides a link between — ; — r— and , , and 

fs{y) " dy 2 dy 



therefore between upward concavity /convexity of concentration profiles and increasing/decreasing of . Since 

js(y) ' - 



d 



i\ j luiiv t;.vi l- y ui Luiiuciiiiaiiun piuinca anu iiii^iccioiiig / ucli ctioing, ui — ' - - 



is always > 0, r— and have the same sign and therefore increasing s allows upward concave concentration 

dy 2 dy fs(y) 

profile, while decreasing £s ^\ allows an upward convex concentration profile. Decreasing profile of £s ^\ in figure (7.b) between 
fs(y) fs{y) 

0.005m and 0.02m allows to predict the near-bed upward convex c(y) profile in figure (8). 

Data of settling velocity w s were derived from measured concentrations for coarse sand and Richardson and Zaki's equation. 
The empirical function f a obtained from these data allows predicting the near-bed upward convex profile. However, the 
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FIG. 8: Predicted time-averaged concentration profile for coarse sand (cj s =6.1 cm/s). Curve, gradient (Fickian) diffusion 
Model with proposed /3(y) and a(y) (= l//s(j/)) functions, x, measurements [131]. 

dimensionless settling velocity function / s decreases from y — 4 cm and at y = 2 cm the decreasing in settling velocity is of 
50%. This seems to be quite larger than would be expected and therefore seems to be outside the range of observed hindered 
settling. 



B. Possible effect of convective sediment entrainment process 



The process of vortex formation and shedding at flow reversal above ripples is a relatively coherent phenomenon. The 
associated convective sediment entrainment process may also be characterized as coherent, instead of a pure disorganized 
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indicated that both convective and diffusive mechanisms are involved in the entrainment processes. In the combined convection- 



diffusion formulation of Nielsen 



3, 



the steady state advection-diffusion equation is given by 



d c 

u) s c + e s (y) uis c F{y) = 

dy 



(34) 



where F{y) is a function describing the probability of a particle reaching height y above the bed. The respective terms in 

1G 



B Possible effect of convective sediment entrainment process 



V CONCLUSION 



Eq. (|34[) represent downward settling, upward diffusion, and upward convection [15( [12| [18| ■ Eq. (|34|l reverts to Eq. (|24l) with 



a(y) = 1 - (co/c) F(y) (35) 
In this equation, a represents the effect of upward convection. The near-bed upward convex profile of figure (8) could be 
predicted by using a parameter ct(y) = l/f s (y) based on (Eq. I30|l . a is given therefore by 



a(y) = l + e*p(-Ui) ( 36 ) 



CONCLUSION 



Field and laboratory measurements of suspended sediments over wave ripples show a contrast between an upward convex 
concentration profiles, time-averaged in semi-log plots, for fine sand and an upward concave profiles for coarse sand. Careful 
examination of experimental data for coarse sand shows a near-bed upward convex profile beneath the main upward concave 
profile. 

The 1-DV gradient diffusion model with a constant settling velocity (w s = cste) and f3 = 1, predicts concentration profile 
for fine sediments but fails for coarse sand. 

The 1-DV gradient diffusion model predicts the main upward concave profile for coarse sediments thanks to a suitable 
/3(i/)-function. The shape of this function was validated by the finite-mixing-length model. 

The model with the resulting sediment diffusivity e 3 and a constant settling velocity is unable to predict the near-bed 
upward convex profile. In order to predict this profile, an additional parameter a is needed. The sediment diffusivity should 
be modulated by a which decreases with y. 

Data of settling velocity w s were derived from measured concentrations of coarse sand with Richardson and Zaki's equation. 
An empirical function for dimensionless settling velocity f s was obtained from these data. The 1-DV gradient diffusion model 
with the proposed function for f s predicts the near-bed upward convex profile. However, the dimensionless settling velocity f a 
decreases from y — 4 cm and at y — 2 cm the decreasing in settling velocity is of 50%. This seems to be quite larger than would 
be expected and therefore seems to be outside the range of observed hindered settling. 

The used parameter a (= 1/ fs) for the near-bed upward convex profile could be related to convective sediment entrainment 
process. Since above ripples the process of vortex formation and shedding at flow reversal is a relatively coherent phenomenon, 
both convective and diffusive mechanisms are involved in the entrainment processes. This seems to be in agreement with the 
combined convection-diffusion formulation of Nielsen 15 ] . 
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V CONCLUSION 



This study was based on experimental data of McFetridge and Nielsen and needs to be tested against other datasets in order 
to allow calibration for predictive purpose. 



References 

[1] Absi R. Comment on 'Turbulent diffusion of momentum and suspended particles: A finite-mixing-length-theory'. Phys 
Fluids 2005;17(7):079101. 

[2] Baldock TE, Tomkins MR, Nielsen P, Hughes MG. Settling velocity of sediments at high concentrations. Coast Eng 
2004;51:91-100. 

[3] Blondeaux P, Vittori G. Vorticity dynamics in an oscillatory flow over a rippled bed, J Fluid Mech 1991;266:257-289. 

[4] Cheng NS. Effect of concentration on settling velocity of sediment particles. J Hydraul Eng, ASCE 1997;123(8): 728-731. 

[5] Conley DC, Falchetti S, Lohmann IP, Brocchini M. The effects of flow stratification by non-cohesive sediment on transport 

in high-energy wave-driven flows. J Fluid Mech 2008;610:43-67. 
[6] Di Felice R. Hydrodynamics of liquid fluidisation. Chcm Eng Sci 2008;50:1213-1245. 

[7] du Toit CG, Sleath JFA. Velocity measurements close to rippled beds in oscillatory flow. J Fluid Mech 1981;112:71-96. 
[8] Dohmen-Janssen CM, Hassan WN, Ribberink JS. Mobile-bed effects in oscillatory sheet flow. J Gcophys Res 
2001;106(C11):27103-27115. 

[9] Fredsoe J. Modelling of non-cohesive sediment transport processes in the marine environment. Coast Eng 1993;21:71-103. 
[10] Graf WH, Cellino M. Suspension flows in open channels: experimental study. J Hydraul Res 2002;40(4):435-447. 
[11] Hsu TW, Jan CD. Calibration of BusingcrArya type of eddy viscosity models parameters. J Waterw Port Coastal Ocean 

Eng, ASCE 1998;124(5):281-284. 
[12] Lee TH, Hanes DM. Comparison of field observations of the vertical distribution of suspended sand and its prediction by 

models, J Geophys Res 1996;101:35613572. 
[13] McFetridge WF, Nielsen P. 1985. Sediment suspension by non-breaking waves over rippled beds. Technical Report No 

UFL/COEL-85/005, Coast Oc Eng Dept, University of Florida. 
[14] Nakato T, Locher FA, Glover JR, Kennedy JF. Wave entrainment of sediment from rippled beds. J Waterway Port Coastal 

Ocean Div 1977;103:83-99. 
[15] Nielsen P. 1992. Coastal bottom boundary layers and sediment transport, World Sci., 324 p. 

[16] Nielsen P, Teakle IAL. Turbulent diffusion of momentum and suspended particles: A finite-mixing-length-theory. Phys 

Fluids 2004; 16(7) :2342-2348. 
[17] Richardson JF, Zaki WN. Sedimentation and fluidisation: part 1. Trans Inst Chem Eng 1954;32:35-53. 



18 



V CONCLUSION 



[18] Thorne PD, Williams JJ, Davies AG. Suspended sediments under waves measured in a large-scale flume facility, J Geophys 
Res 2002;107(C8):3178. 

[19] Thorne PD, Davies AG, Bell PS. Observations and analysis of sediment diffusivity profiles over sandy rippled beds under 

waves, J Geophys Res 2009;114:C02023. 
[20] van der Werf JJ, Ribberink JS, O'Donoghue T, Doucette JS. Modelling and measurement of sand transport processes over 

full-scale ripples in oscillatory flow, Coast Eng 2006;53(8):657-673. 



19 



